

c = 69544098122008365404362084363106703100129828060729951875286300108187409752434091431958070449710849545293894427509731332109603391232669996343824626329922214945852346904863792615793077092573477272770896080330100575733798947163067912513775789041009052964547466568583353539184274864951400026558080635235037201971
c2 =  8943014059616554864780402232625215323717650137702921347505099596041241060991211014413068279308600979202066026580321911322396296877730441259506201930904989
p1 =  9332664074116371272320347910722305705338764655969305966836353509444565128652866496870138156351793824473079367512299950127043084761943406224989567423866071
c3 =  70057129621004770704688083264954773868924488489858841648741006057913061247128697063230635204818097097989773079555422507429472690950271096686354023051578709918278343598395904697714639469429943421928621364157106993940737884253250276769074919155850720804529008723376094183112401317966966358905432385347745890161
p2 =  8345774537268501775456149547604093811537369706242159709058405913660493926215699313366059497170577713800176152866821540944395025269529309431564563983026047
q2 =  10391251387900816431089217882564594263543706933015218000664649428273018708690965953263846227725465978045419390050088292820030415626410428034269853834640767



import libnum
import gmpy2

mp = pow(c2, (p1 + 1) // 4, p1)
print(mp)
print(gmpy2.is_prime(mp))
print(gmpy2.bit_length(mp))

